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Abstract 

A new formulation of chiral fermions on the lattice is presented. It is a version 
of overlap fermions, but built from the computationally efficient staggered 
fermions rather than the previously used Wilson fermions. The construc- 
tion reduces the four quark flavors described by the staggered fermion to 
two quark flavors; this pair can be taken as the up and down quarks in Lat- 
tice QCD. A domain wall formulation giving a truncation of this overlap 
construction is also outlined. 
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1. Introduction 

It is of great interest to find lattice formulations of Quantum Chromody- 
namics (QCD) with exact chiral symmetry. In such a formulation, sponta- 
neous chiral symmetry breaking and the associated spectrum of Goldstone 
bosons (i.e. the light mesons), as well as the emergence of the flavor-singlet 
{rf meson) mass via the index theorem connection between quark zero-modes 
and gauge field topology can be studied directly in the lattice model with- 
out the usual need to extrapolate via chiral perturbation theory. Finding 
such formulations is problematic though 0], and was a long-standing chal- 
lenge in Lattice QCD. An explicit formulation of this kind was finally found 
by Neuberger via the overlap formulation of chiral fermions on the lattice 
[4|. It is a chirally improved version of Wilson fermions, and realizes exact 
chiral symmetry of the Ginsparg- Wilson (GW) type 0,0,0]. However, its 
theoretical attractiveness is offset in practice by the high computational cost 
of implementing it in simulations. 
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In this paper, building on insights from [8j, a new version of overlap 
fermions is presented constructed from the staggered lattice fermion rather 
than the previously used Wilson fermion. The staggered fermion is the most 
computationally efficient lattice fermion formulation known, and this advan- 
tage is expected to be inherited by the new staggered version of overlap 
fermions. The challenge in doing this is that the overlap construction relies 
on specific properties of Wilson fermions that are not shared by staggered 
fermions. This challenge is overcome here using a new theoretical idea of 
interest in its own right: a Wilson-like phase of the staggered fermion theory 
can be created by adding a certain "flavored" mass term. In this new phase, 
the number of quark flavors described by the staggered fermion is reduced 
from four to two, and the staggered fermion acquires Wilson-like properties. 
The latter allow the overlap construction to work as desired. The result- 
ing staggered overlap fermion describes two quark flavors and has an exact 
unflavored GW chiral symmetry originating from the exact flavored chiral 
symmetry of the original staggered fermion. 

However, the flavored vector and chiral symmetries of the quark pair de- 
scribed by the staggered overlap fermion are not exact since the correspond- 
ing symmetries of the original staggered fermion are broken by lattice effects. 
Therefore, exact flavored vector and chiral symmetries only hold for pairs of 
quark flavors when each pair is described by a staggered overlap fermion. 

The two quark flavors described by the staggered overlap fermion can 
be taken as the almost massless u and d quarks in QCD. Then, due to the 
exact unflavored GW chiral symmetry, this description has all the advantages 
of the usual overlap fermions regarding the axial U(l) anomaly and index 
theorem connection between quark zero-modes and gauge field topology. This 
is connected with some of the most subtle and interesting parts of QCD 
physics - not only the rj' mass mentioned already but also the order of the 
QCD phase transition at finite temperature - which are also computationally 
highly demanding to study in lattice QCD (see, e.g., 0]). Finding a new 
lattice approach in which these features are not distorted by lattice effects, 
and which is computationally more efficient than usual overlap fermions, is 
therefore significant. As an indication of the desirability of this, it can be 
noted that high-precision calculation of the rf mass is still an outstanding 
problem in Lattice QCD - the latest state of the art calculation, using Wilson- 
based domain wall fermions (a truncated version of usual overlap fermions), 
has an uncertainly of 15% id . 
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2. Staggered overlap fermion construction 



The staggered lattice fermion field is a one component Grassmann field 
x{x), x = an e aZ 4 , describing 4 degenerate Dirac fermions. The stag- 
gered Dirac operator acting on these fields is D st = r/ M V M where rj^xi 30 ) = 
( — l) ni+ --- +n t j -' 1 x(x) and = ^-(T +fl — T!_ M ), with the parallel transporters 
T±n given in terms of the lattice link variables U^(x) by T +M x(x) = U li (x)x{x+ 
afi) and T_ M = (T +At ) -1 . There is an exact flavored chiral symmetry 

{T 55 ,D st } = (1) 

where r 55 , given by T 55 x(x) = (— l) ni+ --- +n4 x(x), corresponds to 75 <S> 75 in 
the spin-flavor interpretation of staggered fermions [HI]. The first and second 
75 factors act in spinor and flavor space, respectively. Note the properties 

r 55 = 1 and F 55 = r 55- 

Direct application of the overlap construction [3] to staggered fermions 
by simply replacing the kernel operator Dw — M (where Dw is the Wilson- 
Dirac operator) by D st — M does not give anything useful: The exact flavored 
chiral symmetry of the staggered fermion is lost, there is no GW symmetry 
to replace it, and there are no exact zero-modes in general. However, as we 
will now show, the situation changes if, instead of a scalar mass M, we use 
the following "flavor-chiral" mass term: 

M S i = -r 55 r 5 , r>0 (2) 

a 

Here T 5 is the staggered fermion version of the chirality matrix, corresponding 
to 75 <S> 1 up to 0(a 2 ) discretization errors in the spin-flavor interpretation 



111 ]. Consequently, in the spin-flavor interpretation, 

M st ~ 1® 75 + 0(a). (3) 

Explicitly, is given as follows. Use the parallel transporters to define 
= h(T+fj, + T-p), and let C = {C\C2CzCi) syrn denote the symmetrized 
product of Ci, C2, C3, C4. Then T 5 = r] 5 C where 775 = 771772^3^4- (Explicitly, 
V5X( X ) = (~ l) ni+n3 x( x )-) Note that T 5 is hermitian and commutes with r 55 ; 
therefore M st is hermitian. The effect of this mass term in 

D st ~ M st (4) 

is to split the low-lying modes into branches with approximately definite pos- 
itive and negative flavor-chirality, giving them masses — 1 and -, respectively. 
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Figure 1: Spectrum of D st — M st in the free field case 12[. 



This is shown in the free field case in Fig. 1, taken from [12j. In this case, 
there are 16 zero-modes of D st with momenta pa = -A, A^ e {0,1}. On 
the vectorspace V spanned by these modes we have (r 55 r 5 ) 2 = 1, so there 
is a decomposition V = V + © V- with r 55 r 5 = ±1 on V±, and it can be 
shown that dim(V±) = 8. Hence D st — M st = =F^ on V±; this shows how the 
real eigenvalues arise in the free field case. Moreover, the operators f]fj,T +fl 
and rj^T-n making up D st coincide on V and can be shown to commute with 



on V. Hence their representations T M on V (given in [TlJ by the 16 x 16 
matrices (T^)ab), which form a representation of the Euclidean Dirac alge- 
bra (llj . decompose into representations on V±. From this it can be seen 
that, from a low momentum viewpoint, the 4 Dirac fermion species described 
by the staggered fermion are split into two pairs of Dirac fermions, with the 
pair with positive (negative) flavor-chirality having mass — - (+-). 
Upon inserting the kernel D st — M st into the overlap formula, 



-(l + (D s 

a \ 



D so = - ( 1 + (D st — M st 



^(D st -M st y(D st -M st ). 



(5) 



the negative mass modes are converted into the physical modes of a massless 
fermion, while the positive mass modes are converted into modes with heavy 
masses ~ 1/a which decouple in the continuum limit, just as in the Wilson 
case. Thus only two of the original four quark flavors survive in the resulting 
overlap fermion. Furthermore, since the physical modes have approximately 



4 



definite positive flavor-chirality, T 5 = r 55 on these modes up to 0(a 2 ) effects. 

The staggered overlap Dirac operator fl^D satisfies the following GW re- 
lation: 

{r 55 , D so } = ^D so r 55 D so (6) 

Note that this is now an unflavored chiral symmetry, since, as remarked, 
= r 5 +0(a 2 ) on the physical modes of the staggered overlap fermion. The 
GW relation can be readily verified after noting that the staggered overlap 
operator can be expressed as 



Dso = - i + r 55 ^= (7) 




with the hermitian operator H st given by 

H st = r 55 {D st - M st ) = r 55 D st - ~r 5 (8) 

(X 

It also follows from ([7]) that D\ Q = r 55 .D so r55, which implies that the non-real 
eigenvalues of D so come in complex conjugate pairs. Since the only possible 
real eigenvalues (which must lie on the GW circle) are and 2r/a, it follows 
that the fermion determinant det(D so ) is real and non- negative, as required 
for lattice QCD simulations. 

A required property of the lattice Dirac operator is locality. It holds 
for D so ; this can be proved analytically by the same argument as in the 



Wilson overlap case [141 ]. using the staggered bound mentioned in [8], when 



the lattice gauge fields satisfy an admissibility condition; the details are in 



13l | . Locality of D so in realistic gauge field backgrounds has been shown 



numerically in [12 



Computational efficiency of D so was also studied in [12j : preliminary 
numerical evidence indicates a speedup of order 2-3 compared to the Wilson- 
based overlap Dirac operator in realistic gauge field backgrounds. 



3. Pairs of exact zero-modes and index theorem 

We now verify that D so has pairs of exact chiral zero modes with index 
determined by gauge field topology in accordance with the 2-flavor index 
theorem when the lattice gauge field background is not too rough. Due to 
the GW relation fl^]), the index is well-defined as index (D so ) = n + — n_ 
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where n± is the number of independent zero-modes of D so with ± chirality 
under r 55 (which is the same as T 5 on these modes up to 0(a 2 ) effects). A 
general formula j^j gives 



mdex(D so ) = -~Tr . (9) 




We connect this to the hermitian staggered operator studied previously in 
[ijj], denoted by H st below: 

Theorem. For every value of m, the hermitian staggered operators 
H st (m) = T 55 D st - mT 5 , H st (m) = iD st - mT 5 

have the same eigenvalue spectrum. 

The proof of the theorem is deferred to 13|. It can be understood intu- 
itively from the fact that H st (m) arises from H st (m) by a change of repre- 
sentation of the staggered sign factors — > iT^rj^ in D st . 

In light of the theorem we can replace H st in by H st (m) with m = -; 
here r is the parameter in (TJJ). But then the results of [s| give 

mdex(D so ) = §mdex(D st ) = 2Q (10) 

in sufficiently smooth gauge field backgrounds of topological charge Q and 
with - lying in a suitable range. The last equality in ( FlOl) was confirmed in 
the numerical study in 

4. Stability of the massless phase under radiative corrections 

It is important to consider the symmetries of the staggered overlap fermion 
since these need to be sufficient for renormalizability and stability of the 
massless 2-flavor phase under radiative corrections. Besides the replacement 
of the flavored chiral symmetry (pQ) b y t he GW symmetry ([6]), the usual 
staggered fermion symmetries listed in all hold for the staggered overlap 
fermion, except for axis reversals and shift transformations. Under the latter 
xr 55 r 5 x — > — x^I^X) and therefore \D S oX is invariant when these trans- 
formations are combined with the parameter flip r — > — r in (j2]). The same 
invariances must hold for all counterterms that arise. Hence the counterterms 
must be invariant under all the usual staggered fermion symmetries except 
([Q) , and except for a possible sign change under axis reversals and shift trans- 



formations. It is easy to show [13| that the most general local mass-dimension 
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4 terms with these properties are (up to proportionality, and modulo 0(a) 
terms )0 

1 _ 1 _ 

xDstx , x{r 5 5r 5 , D st }x , -x^55^5X , -xx ( n ) 

a a 

Since _D so is local, all counterterms are also local, and the theory is renor- 



malizable, provided the power-counting theorem holds [15|. It holds in the 



Wilson overlap case [16[ and is plausible in the present case. Assuming this, 
stability can be established to all orders in the loop expansion as follows. 
The GW relation (jBJ) leads to 

r a (D-\q)) + (D-\ q ))r BB = ±r 8li (12) 

where (D~}(q)) = (D so (q) + S(g)) _1 with E(g) being the fermion self-energy. 
Here and in the following, X(q) denotes the free field momentum repre- 
sentation of the staggered operator X; it is a linear map (16 x 16 ma- 
trix) X(q) : V — > V on the zero- mode vector space V defined earlier, with 
g M G [— j-, ^]. The momentum representation of r 55 is a constant matrix, 
and r 55 r 5 (g) = ±C(aq)l on V± where C(aq) = Ylij, cos ( a( lti.) = 1 + 0(a 2 q 2 ). 
In light of the preceding, at one loop £(?) has the form 

CiD st (q) + c 2 T 55 T 5 (q)D st (q) + -r 55 r 5 (g) + ^ (13) 

a a 

up to irrelevant terms, with each Cj = Cj(a 2 q 2 ) diverging no worse than 
logarithmically for a — > 0. Note that all terms in ( fTBl map V± to itself, 
as does D so (q). From (jT2 |) -( |T3l) and the fact that D so (0) = on V + it is 
straightforward to show S(0) = on V+. (We omit the details; an analogous 
result was found in the Wilson overlap case in [l6[ [13]). This implies c 4 = — c 3 
in ( TT3|) . so on V + the last two terms in ( TT3|) combine to give an irrelevant 
term ^-(C(aq) — 1) which vanishes ~ alog(a) for a — > 0. The first two 
terms in (TT3T) act on V± as (ci ± c 2 C(aq , ))(r^ sin(ag M )). It follows that, for 
a -)• 0, ,D so (g) + E{q) = T+iq^l + c^a 2 /) + c 2 {a 2 q 2 )) on V+. The self- 
energy is hereby seen to have the usual effect of a logarithmically divergent 
wave function renormalization on the massless physical modes. Then, after 
renormalizing the theory at one loop, the argument above can be repeated 
to get the same result at two loops, and so on. In this way stability is seen 
to hold to all orders in the renormalized theory. 



1 This is with the lattice gauge fields included. In the free field case r 55 r 5 and D st 
commute and the second term in (|11J simplifies. 
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5. Wilson-like nature of the staggered overlap fermion kernel 

A deeper understanding of the staggered overlap construction can be 
obtained by clarifying the Wilson-like nature of its staggered fermion kernel. 
To this end we introduce the following 2-flavor Wilson-like version of the 
staggered Dirac operator: 

D sW = D st + W st , W st = -(l-V b5 V 5 ). (14) 

a 

The term W st plays an analogous role to the Wilson term in the Wilson- 
Dirac operator: it decouples the negative flavor-chirality modes by giving 
them mass 2r/a while keeping the two positive flavor-chirality modes as 
the physical modes. Hence D s w describes two physical quark flavors on 
which Y§ = up to 0(a) effects. In light of this we can obtain a 2- 
flavor overlap fermion by taking D s w — M with M = rp/a, p G (0,2) as 
the kernel in the usual overlap construction. For p — 1 this is precisely the 
2-flavor staggered overlap Dirac operator D so we constructed above in (JSJ), 
since D s w — - = D st — M st . But now we see that it can be generalized to any 
p G (0, 2). Furthermore, the role of the parameter r in the staggered overlap 
construction is hereby clarified: it is analogous to the Wilson parameter in 
the usual overlap construction based on Wilson fermions. This is surprising, 
since r/a in ([2]) initially appears to be analogous to the mass parameter 
M = p/a in the Wilson case. 

Starting from Wilson fermions, a domain wall fermion formulation can be 
constructed [18] which gives a truncation of the overlap fermion construction 



19| . The same can now be done with staggered fermions simply by replacing 

D w -»■ D sW , 75 r 55 (15) 

in the previous Wilson-based constructions. The lattice Dirac operator for 
the staggered domain wall fermion in 5 dimensions is then 

D sdw = D sW -M + T 55 d s , M = ^ , p E (0, 2) (16) 

where s G [0, L\ is the lattice coordinate of the 5th dimension and d s = 
P + d { s +) + P_dt ] with P± = |(1 ± r 55 ) and d { s +) , <9 S H being the forward 
and backward finite difference operators^ Boundary conditions on x(x, s) 



2 This is completely unrelated to a previous proposal for staggered domain wall fermions 
in 120| which describes 4 rather than 2 flavors. 
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at s = and s = L are introduced analogously to the Wilson-based case; 
then the staggered domain wall construction is seen to be a truncation of the 



staggered overlap construction by essentially the same argument as in [19 



The details will be given in [13 



6. Concluding discussion 

Underlying the staggered overlap and staggered domain wall fermion con- 
structions in this paper is a staggered version of Wilson fermions obtained 
by adding a flavored mass term to the staggered fermion action. The idea 
of considering staggered fermion with flavored mass term is not new in it- 



self; this was already done many years ago by Golterman and Smit in [11 
Rather, the new idea here is that, for a certain choice of flavored mass term, 
namely the Wilson term W st in ( |14p determined by the "flavor-chiral" mass 
term M st in (j2J), one can use the staggered operator r 55 for the role of 75 
since it coincides with 75 <g> 1 on the 2 physical flavors up to 0(a) effects. The 
significance of this is that the properties 75 = 1 and 75— hermiticity, which 
are crucial in the usual overlap and domain wall constructions, continue to 
hold in the staggered setting with 75 — > T 55 . This would not be true if one 
tried to use the direct staggered analogue T 5 of 75. 

The flavor-chiral mass term used in this work, which reduces the 4 stag- 
gered fermion flavors to 2 flavors, is not the only possibility. Another flavor- 
chiral mass term which reduces the number of flavors to 1 was subsequently 
proposed in [2l| after the first version of this paper appeared on the Arxiv. 
However, it breaks more of the staggered fermion symmetries than the present 
proposal, namely some of the lattice rotation symmetries. The consequences 
of this for radiatively generated counter-terms and fine-tuning requirements 
remains to be seen. The 1-flavor version of staggered overlap fermions has the 
advantage that all the flavored vector and GW-chiral symmetries are exact. 
On the other hand, for problems where only the unflavored chiral symmetry 
is important, e.g. for computing the 7/ mass, the present 2-flavor formulation 
has advantages: it preserves lattice rotation symmetry and gives "two quarks 
(the u and d) for the price of one" in lattice QCD simulations. 
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